Measurement of P M £ in Polarized Muon Decay 
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The quantity P££, where £ is one of the muon decay parameters and P* is the degree of muon po- 
larization in pion decay, has been measured. The value — 1.0003T0.0006 (stat.) ±0.0038 (syst.) 
was obtained. This result agrees with previous measurements but is over a factor of two more pre- 
cise. It also agrees with the Standard Model prediction for P£ £ and thus leads to restrictions on 
left-right symmetric models. 

PACS numbers: 13.35.Bv, 14.60.Ef, 12.60.Cn 



I. INTRODUCTION 

In the Standard Model (SM) of particle physics, pos- 
itive muons decay via the weak (V — A) interaction into 
positrons plus neutrinos: p^evv through a virtual state 
involving W vector bosons. More generally, the ampli- 
tude for muon decay can be described in terms of a local 
decay matrix element, which is invariant under Lorentz 
transformations: 



M= ^f E ^<ee|r>)(77|r> M >, (i) 
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where the g]^ specify the scalar, vector, and tensor cou- 
lings between /i-handed muons and e-handed positrons 
1]. In the SM g\ h = 1, and all other coupling constants 

are zero. 

The differential decay spectrum Q of the e + emitted 
in the decay of a polarized fi + can be described by four 
parameters - p, S, 7/ and £ - commonly referred to as 
the Michel parameters, which are bilinear combinations 
of the coupling constants. In the limit where the positron 
and neutrino masses are neglected, and radiative correc- 
tions |jj are not explicitly included, this spectrum is given 
by: 



where 9 is the angle between the muon polarization and 
the outgoing positron direction, x — E e /E max , E max — 
(to 2 + to 2 )/2to m = 52.83 MeV, and P M is the degree of 
muon polarization. The fourth parameter, ry, appears in 
the isotropic term when the positron mass is included 
in the analysis. In the SM, the Michel parameters take 
on the precise values p = S = 0.75, £ = 1, and 77 = 0. 
The parameter £ expresses the level of parity violation in 
muon decay, while S parametrizes its momentum depen- 
dence. 

In this experiment P M is the magnitude of the p, + po- 
larization along the beam axis at the time of muon decay. 
Surface p + [J], which are muons produced from 7r + de- 
cays at rest, have a polarization of magnitude P£, with 
a direction antiparallel to their momentum, given by a 
generalization of Eq. 2] for scmileptonic decays. In the 
SM with massless neutrinos P£ = 1. In this experiment 
P M £ is determined from the positron spectrum, and then 
P^£ is obtained using the measurements of the muon 
trajectories. 

SM extensions involving right-handed interactions [f| 
require deviations from pure V—A coupling that can alter 
P M £. Four probabilities Q e/i (e, p, = R, L) for the decay of 
a /j-handed muon into an e-handed positron are given by: 
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II. EXPERIMENTAL PROCEDURES 
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where S efi = 1 for e = p and <5 £A1 = for e =^ /i. The 
probability: 
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sets a model independent limit on any muon right-handed 
couplings El • A recent review of muon decay is pre- 
sented in 0| . 

A precision measurement of muon decay can place lim- 
its on left-right symmetric (LRS) models The LRS 
models contain four charged gauge bosons (W^, W^), 
the photon, and two additional massive neutral gauge 
bosons. The W\ and Wi masses are m\ and mi respec- 
tively, and the fields Wl and Wr are related to the mass 
eigenstates W\ and Wi through a mixing angle £■ In 
these models both V— A and V+A couplings are present, 
and parity violation appears because of the difference in 
the mass of the vector bosons. In the SM the weak in- 
teraction contains only two left handed vector bosons, 

wi- 

In the general LRS model (5j, 



£ ~ 1 - 2 
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9l mi 
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(5) 



where and g^ are the right- and left-handed gauge 
couplings. The manifest left-right symmetric model 
makes the additional assumptions that gn = g^ and 
that the left- and right-handed quark mixing matrices 
are identical. In this case, P™ can also be expressed in 
terms of mi /mi and £, and one obtains 



mi 
m 2 
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mi 
m 2 



(6) 



Prior to TWIST, the most precise direct measurement 
of P*£ was 1.0027 ± 0.0079 (stat.) ±0.0028 (syst.) @, 
in agreement with the SM. A similar value has been 
measured using muons from kaon decay Recently 
the TWIST collaboration reported new measurements 
of p = 0.75080 ± 0.00032 (stat.) ± 0.00097 (syst.) ± 
0.000023 (r?) [Sand 5 = 0.74964 ± 0.00066 (stat.) ± 
0.00112 (syst.)|n|- Using the result P^8/p > 0.99682, 
at the 90% confidence level 12], along with the TWIST 
measurements of p and 5, an indirect limit on was 
determined to be 0.9960 < P^ < 1.0040 (90% con- 
fidence level) |lll |. In this paper a new measurement of 
is reported. 



In the present experiment, highly polarized surface 
muons were delivered, in vacuum, to the TWIST spec- 
trometer by the M13 channel at TRIUMF 0. The 
properties of the surface muon beam are a typical rate 
of 2.5 kHz, a momentum of 29.6 MeV/c, and a momen- 
tum bite, Ap/p w 1.0% FWHM. Beam positrons with a 
typical rate of 22 kHz, and the same momentum as the 
muons are also a feature of the surface muon beam. 

A low pressure (8 kPa dimethyl ether gas) removable 
beam monitoring chamber [lfij provides information on 
the muon beam before it traverses the fringe field of the 
solenoid. The chamber consists of two 4 cm long modules, 
one to measure the position and direction of the muon 
beam in the horizontal (x) direction, and the other for 
measurements in the vertical (y) direction. The location 
of the final quadrupole of the M13 channel, the beam 
monitor, muon ranging gas degrader, trigger scintillator 
and the TWIST detector are shown in Fig. ^ The beam 
monitor is inserted for measurement of the beam proper- 
ties, and removed during data collection. 
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FIG. 1: Location of the last beamline quadrupole, beam 
monitor, gas degrader, trigger scintillator, and the TWIST 
solenoid. 



Changing the angle of the muon beam relative to the 
magnetic field axis gives rise to a change in the polar- 
ization. The solenoid field is found to interact with the 
iron of the beamline magnets, such that the muon beam 
is deflected off axis. The available Ml 3 channel magnets 
could only partially alleviate this deflection. Figure Ufa) 



3 



shows that the beam is centered in x, but not in y. The 
measured muon beam distribution in position and angle 
at the beam monitor, as shown in Fig. |2 is input into 
a simulation to calculate the average depolarization of 
the muons from the location of the beam monitor to the 
high-purity Al muon stopping target at the center of the 
TWIST detector (Fig. gj. The angular distribution in x 
is shown in the figure, but similar distributions are mea- 
sured in y and these are also included in the simulation 
(lij. The RMS size of the beam is 0.6 cm in both x and 
y, and the average angular spread of the beam is 10 mrad 
in x and 12 mrad in y. 
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FIG. 2: (colour online) Beam measurements projected to the 
center of the beam monitor at each 1 mm by 1 mm bin in y 
versus x. (a) Muon beam intensity profile, (b) Mean angle in 
x (6 X ). 

The muon beam is transported in vacuum until it is 
inside the solenoid field, where it then passes through the 
gas degrader and trigger scintillator. The gas degrader is 
a 21.67 cm long (along z) chamber containing a mixture 
of He and CO2 that can be adjusted to stop the muons 
in the stopping target. The plastic (Bicron BC408) disk 
shaped trigger scintillator has a thickness of 195 /im, a 
radius of 3 cm, and is located 80 cm upstream of the 
central stopping target. 



The materials that the muons pass through, starting 
from the production target to the center of the stopping 
target, are summarized in Table [I] 



Material 


Thickness (mg/cm 2 ) 


j_yegiaoei alio vaeLiLim 1011a 


1 1 Q1 


He/C02 degrader 


1.95 to 42.80 


Air gap (2.82 cm) 


3.65 


Muon scintillator 


20.12 


Scintillator wrap 


3.03 


Cradle window 


0.88 


PC module 


9.46 


Dense stack 


13.41 


Seven UV modules 


27.80 


He/Na (63.8 cm) 


12.25 


Half target module before target 


4.31 


Half of 71 /im Al target 


9.59 


Total to center of Al target 


118.36 to 159.21 



TABLE I: Estimates of the material thicknesses in mg/cm 2 
that the muon penetrates from the production target to the 
center of the muon stopping target. Surface muon momenta 
are such that this approximately matches the muon range for 
a degrader containing 50% CO2. 



The TWIST spectrometer |l3j is designed to measure a 
broad range of the positron spectrum from muon decays 
in a stopping target, allowing the simultaneous extrac- 
tion of the spectrum shape parameters. The spectrom- 
eter consists of 12 very thin high-precision proportional 
chamber (PC) planes and 44 drift chamber (DC) planes, 
perpendicular to the axis of a solenoid producing a mag- 
netic field of 2 T. The upstream half of the spectrometer 
(22 DCs and 6 PCs) is shown in Fig. [31 and the down- 
stream half of the detector is mirror symmetric with the 
upstream. 

The gas particle detectors used by TWIST are 
mounted in a leak tight aluminum cradle that is filled 
with 3% N2 and 97% He gas. Helium, rather than air, 
is used to reduce the scattering between measurement 
points and to allow the low energy muons to reach the 
stopping target at the center of TWIST. The spacing 
of the PCs and DCs is maintained by precision ground 
spacers made of a Russian engineered material Sital CO- 
113M which has a very small thermal coefficient of expan- 
sion. The spacers are pushed together by four pneumatic 
cylinders with a force of 1470 N each to ensure they are 
touching and do not move. 

The DCs are used to obtain a precise measurement of 
the position of the e + as it traverses the chambers. To 
accomplish this a slow drift gas, dimethyl ether (DME), 
which has a small angle between the drift field and 
the electron drift direction in a non-zero magnetic field 
(Lorentz angle), is used. Each DC plane consists of 80 
gold-plated tungsten anode sense wires of 15 /im diame- 
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FIG. 3: (colour online) Side view of the upstream half of the 
TWIST spectrometer planar chambers and support structure. 

ter, spaced by 4 mm ^||. The cathode foil walls are 6.35 
/an of aluminized Mylar, nominally placed 2 mm from 
the wires. The drift time for each signal is used along 
with a space-to-time relationship (STR) that is calcu- 
lated by a GARFIELD simulation |l7j to improve the 
determination of the momentum and angle of the decay 
positrons. 

The purpose of the PCs is to have reasonably fast tim- 
ing information to help with pattern recognition. The 
width of the time to digital converter (TDC) signal from 
the PCs is used to discriminate muons from positrons. 
The PCs are similar to the DCs, with the exception of 
the following: the PCs use CF^isobutane gas, a wire 
spacing of 2mm, and 160 wires in each plane. The PC 
modules consist of four wire planes. One PC module is 
placed at either end of the detector stack, and the spe- 
cialized target module in the center of the detector stack 
contains two PCs on either side of the muon stopping 
target. 

Collection of an event is triggered by the muon scin- 
tillator whose threshold is set high enough to reduce the 
beam positron trigger rate but low enough to trigger on 
most beam muons. The fraction of triggers due to beam 
positrons is only 10%. Positrons in the beam are easily 
removed from the data since they leave signals through 
the full length of the detector. 

Events in TWIST are recorded by LeCroy Model 1877 
TDCs. The trigger and TDC read-out are set up to 
record signals in 500 ps time bins from 6 /zs before to 10 
/Lts after a muon passes through the trigger scintillator. 
For each wire signal, the TDC records a pulse start and 
stop, and thus can be used to calculate a pulse width. In 
this configuration up to eight wire signals can be recorded 



for each wire in any triggered event |l8| . A fixed blank- 
ing time of 80 jj,s was imposed after each accepted trigger. 
This blank is sufficiently long to allow each TDC to finish 
conversion. To insure that all TDCs have fully converted 
their hits, each one was required to produce a header 
specifying how many hits were converted even when zero 
hits were received. 

For this measurement a 71 ±1 fim thick, 99.999% pure, 
Al target is used to stop muons. The tracks of the se- 
lected muons are required to have at least one signal in 
the PC immediately before the Al target and no signals 
in the PC immediately after the Al target. The average 
z from the last chamber plane fired by the muons (z^) is 
used in a feedback loop to control the fractions of He and 
CO2 in a gas degrader, in order to stop 97.0 ± 1.5% of 
the selected muons in the Al target at the center of the 
solenoid. Some muons that stop in the CF^isobutane or 
Mylar foils in the vicinity of the Al target could not be 
separated and wer included in the measurement. 

III. DATA ANALYSIS 

TWIST determines the Michel parameters by fitting 
two-dimensional distributions of reconstructed experi- 
mental decay positron momenta and an gles with distri- 
butions of reconstructed simulated data |ll| . 

The nominal fiducial region adopted for this analysis 
requires p < 50 MeV/c, \p z \ > 13.7 MeV/c, p T < 38.5 
MeV/c, and 0.50 < |cos0| < 0.84. The fiducial cuts, 
while intentionally chosen to be conservative, are related 
to physical limitations of the TWIST detector. The 50 
MeV/c momentum cut rejects events that are near the 
region utilized in the energy calibration. It is also im- 
portant to avoid the region very close to the end point 
to minimize the sensitivity of the Michel parameter fits 
to details of the simulation that may affect the momen- 
tum resolution. The longitudinal momentum constraint 
eliminates events with helix pitch near the 12.4 cm pe- 
riodicity in the wire chamber spacing. The transverse 
momentum constraint ensures that all decays are well 
confined within the wire chamber volume. The angu- 
lar constraint removes events at large |cos#| that have 
worse resolution and events at small |cosf?| that expe- 
rience large energy loss and multiple scattering. These 
limits were fixed early in the analysis. The value of P££ 
was found to change by less than 0.0001 when the fidu- 
cial boundaries were moved by ±2% in momentum cut 
values and ±10% in | cos6>| cut values. 

The decay positrons spiral through the chambers pro- 
ducing signals on the wires, which are recorded by the 
TDCs. These helical tracks are subsequently recon- 
structed and analyzed to determine the positron energy 
and angular distributions. Determining the positron mo- 
mentum and angle is done with a x-squared fit to a helical 
track that includes the drift time information from each 
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cell. The helix fits also include positron multiple scatter- 
ing in the y 2 calculation using the procedure from |l9| . 
The efficiency of the track fitting is >99.5% within the 
nominal fiducial region used for spectrum fitting [Tol. Hlf . 
The number of degrees of freedom (DOF) varies with 
track angle, since higher angle tracks will pass through 
more DC cells per plane. For the average DOF, the most 
probable x 2 /DOF is 29.7/28 = 1.06, but there is a long 
tail to larger values because the non-Gaussian tails of the 
multiple scattering distributions are not treated prop- 
erly in the \ 2 calculation and because we only use ap- 
proximate STRs. The momentum resolution is typically 
100 keV/c, and the cos# resolution is about 0.005 [Io| . 

In event selection, the PC and DC signals are exam- 
ined to identify events in which the muon stopped in the 
target, then decayed at least 1.05 (is, and no more than 
9 fis, later. The delay insures that the PC and DC sig- 
nals associated with the muon and decay positron do not 
overlap. Events are rejected if a second muon enters the 
spectrometer, or if a beam positron passes through the 
spectrometer within 1.05 fis of either the muon arrival or 
decay time. Additional cuts include the muon flight time 
through the M13 beam line and a requirement that the 
muon stopping location be within 2.5 cm of the detec- 
tor axis. All events that pass these cuts are analyzed to 
reconstruct the decay positron kinematics. 

After track fitting, ~ 2% of the events contain addi- 
tional tracks in coincidence with the decay. Extra tracks 
can arise from beam particles that are not resolved in 
time, events that scatter within the detector leading to 
two reconstructed track segments, and events that in- 
clude delta rays or decay positrons that backscatter from 
material outside the detector volume. Algorithms were 
developed to identify events with backscatters and reject 
coincident beam particles [ll|. 

The energy calibration of the decay positrons is ob- 
tained from a fit to the endpoint of the spectrum. The 
fit for the endpoint is done separately from the fit for 
the muon decay parameters, and to avoid bias the data 
from the region of momentum used in the endpoint fit 
are excluded from the determination of the muon decay 
parameters. The endpoint fit function is a slope with an 
edge (/(x) = ax + b for x < 1, and f(x) = for x > 1) 
convoluted with a Gaussian a. The parameters a, 6, and 
a all depend on cos 8 [2Cj . The end point of the muon 
decay spectrum and sections of the 2-dimensional end 
point fit function for the bins within the fiducial region 
with the smallest upstream and downstream angles are 
shown in Fig. 01 The difference in yield between up- 
stream and downstream emphasizes the asymmetry of 
polarized muon decay. The corrected momentum p oc is 
given by: 



Pcdg 



cost 



(7) 



reconstructed cosine of the decay positron angle, p e dge is 
the maximum positron momentum, (3 defines the momen- 
tum scale related to the magnitude of the spectrometer 
magnetic field, and a = {a u , ad) is the zero angle energy 
loss for upstream (u) or downstream (d) decay positron 
tracks. This simple form was chosen because of the pla- 
nar geometry of the wire chambers, such that the amount 
of material the positron passes through increases linearly 
as 1/ cos9. 



-0.84 < cosB < -0.82 




51.5 



52 52.5 53 

Momentum (MeV/c) 



53.5 



where p rcc is the reconstructed momentum, cos 6 is the 



FIG. 4: (colour online) Sections of the 2-dimensional end 
point fit function of the muon decay spectrum for the bins 
within the fiducial region containing the smallest upstream 
and downstream angles. The data are shown as a solid-line, 
the matching simulation set is shown as a dashed-line, and 
the endpoint fit functions are shown as smooth curves. 

The endpoint fit parameter (3 is highly correlated with 
c*sum = <^u+ctd- In the simulation the momentum is mea- 
sured without bias at the few keV/c level, the magnetic 
field is measured with an NMR probe, and the magnetic 
field map has been determined to better than 0.2 mT 
in the tracking region. For this reason the endpoint fits 
are done with the value of (3 set to zero and assigned an 
uncertainty consistent with the momentum fit and field 
map accuracy. A 12 keV/c difference between the data 
and simulation cv SU m values was observed, and this is cor- 
rected by applying the energy calibration. 

The TWIST simulation is based on GEANT 3.21 
The simulation contains virtually all the components of 
the spectrometer with which a muon or a decay positron 
could interact. The output exactly mimics the binary 
files generated by the data acquisition system. Details of 
the simulation have been presented previously 0, E3 • 

In bringing the muon beam to a stop in the Al target 
at the center of the TWIST spectrometer, the muons are 
depolarized by the combined effect of multiple scattering, 
interaction with the fringe field of the spectrometer, and 
interactions when stopped in the high purity Al target. 
Thus, the polarization of the muon with respect to the 
z axis when it decays is lower than its polarization with 
respect to the muon momentum when it was produced in 
pion decay. To obtain an absolute measurement of P££,, 
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data are fit to a simulation that includes effects of fringe 
field depolarization and material depolarization. 

The main factors that influence the difference between 
P£ and Pfj, are the transport of the muon spins in the 
various regions of magnetic field and models for the 
muon depolarization in materials. Transport of the muon 
spins from the beam measurement location to the muon 
stopping target is done using a classical fourth order 
Runge-Kutta using the Nystrom algorithm [2^], of the 
Bargmann-Michel-Telegdi |2^| equation. The inherent 
accuracy of this numerical integration is excellent, but 
depends on the accuracy on the knowledge of the in- 
put beam parameters and the magnetic field map. For 
this reason, the beam monitoring chamber measurements 
were used to generate the muon beam for the simulation. 

Ideally the magnetic field is meant to be uniform and 
aligned with the axis of the muon beam momentum (and 
anti- aligned with the muon spin), which is labeled as the 
+z coordinate in TWIST. The field shaping elements 
of the magnet were input into a finite element analy- 
sis which models the magnetic field in order to create 
a field map that included the radial components of the 
field as well as the longitudinal components at each of 
the positions at which the longitudinal component was 
measured. 

The magnetic field of the solenoid was mapped using 
a specially constructed mapping tool which used Hall 
probes, and NMR magnetometers to measure the mag- 
netic field. The component of the magnetic field along z 
was mapped in the uniform field region and in the beam 
monitoring entrance region. A field map was generated 
that matches the measurements to 0.2 mT in 2.0 T in the 
uniform region, but deviates up to 4% in the entrance re- 
gion at z < —200 cm. The effects of the steel from the last 
two quadrupoles of the MI3 beamline were studied using 
field maps from the finite clement analysis that included 
or excluded the iron of the last three beamline elements. 
The map produced from this study was used to estimate 
a systematic error on P^t; due to the uncertainty in the 
magnetic field in the fringe field region. 

A blind analysis was implemented by utilizing hidden 
Michel parameters ph,5h, and to generate the simu- 
lated decay rate. The decay rate can be written as: 



d 2 T 



dxd(cos ( 



d_ 



d 2 T 



dxd(cos ( 



AA, 



because the decay spectrum is linear in the shape pa- 
rameters. The simulation spectrum was fit to the data 
spectrum by adjusting A/?, A£, and A£<5. 

An alternate analysis scheme, used only to compare 
relative polarizations, was developed using an integral 
asymmetry defined as the difference between the number 
of forward and the number of backward decays divided 
by their sum. To obtain a polarization estimate the for- 
ward and backward sums were done inside the fiducial 



region described earlier in this section and normalized 
using integrals of Eq. (2Jwith the SM values of the Michel 
parameters inserted. 



IV. EVALUATION OF SYSTEM ATICS 

The leading systematic uncertainties in this measure- 
ment of arise from the potential sources of muon de- 
polarization. These include depolarization due to multi- 
ple scattering in the production target and a 3 /jm beam- 
line foil, fringe field depolarization, and interactions with 
material while the muon is propagating through the de- 
tector and after stopping. 

The depolarization in the production target is due to 
multiple scattering of the muons while exiting the target. 
An estimate of the depolarization in the small angle ap- 
proximation is (6»^ c s c ) 2 /2, where 6^ c s e is the width of the 
multiple scattering angular distribution |(|. The muons 
in the beam arise from a maximum depth of 0.003 cm of 
graphite, which produces a depolarization of 0.2 x 10~ 3 
which is taken as the systematic uncertainty in due 
to multiple scattering in the production target. 

The systematic uncertainty in the fringe field depolar- 
ization is estimated from the different settings used in 
data taking for the second dipolc element (B2) in the 
M13 channel. The beam parameters were measured for 
two different B2 settings, both before and after the data 
collection. The relative changes in angle and position 
between the nominal B2 value (94.4 mT) and B2+0.5% 
settings are similar for the two periods, but the absolute 
numbers for the average beam angles are quite differ- 
ent. This could be due to changes in the performance of 
the beam monitoring chamber or to its alignment to the 
beamline. To determine the sensitivity of the polariza- 
tion to beam position and angle, a simulated beam was 
scanned in position and angle and the polarization was 
found to depend quadratically on the input variables. Us- 
ing this parameterization, the predicted polarizations for 
the four characterization runs are shown in Table ITT1 The 
beam for the third row of Table [H] is illustrated in Fig. 
12 The larger of the differences in predicted polarization 
for a given B2 setting (0.0033) is adopted as an estimate 
of the uncertainty due to limits of reproducibility. 

Uncertainties due to deconvolution of the beam an- 
gle measurement, modeling of the shape of the solenoid 
fringe field, and beam size reproduction also contribute 
to the final quoted systematic uncertainty of 0.0034 due 
to fringe field depolarization. 

The depolarization of the muons while they propagate 
through the detector and interact with the detector ma- 
terials is negligibly small for muons with kinetic energy 
between 4 MeV and 100 keV ■ Simulations show that 
greater than 99% of muons stopping in our Al target 
have more than 100 keV kinetic energy. Muonium for- 
mation is thus suppressed by ensuring that the majority 



7 



TABLE II: Average beam positions and angles from beam 
monitoring measurements taken at different times, along with 
the simulation estimates of the muon polarization. The first 
(last) two entries are from before (after) the data collection. 



half the difference between the Gaussian and exponential 
extrapolations. 



B2 
(mT) 


X 

(cm) 


0.c 

(mrad) 


y 

(cm) 


9y 

(mrad) 


psim 


94.4 


0.07 


-5.9 


0.97 


7.0 


0.9929 


94.9 


0.85 


-1.1 


0.87 


-5.0 


0.9955 


94.4 


0.06 


-6.7 


0.73 


-11.2 


0.9941 


94.9 


0.94 


-1.5 


0.64 


-19.2 


0.9922 



of the muons have sufficient energy entering the Al tar- 
get. For those muons that do not stop in the Al target, 
the Paschen-Back effect minimizes but does not eliminate 
the depolarization due to muonium formation. Most of 
the muons stop in the high-purity Al target, where they 
can interact with conduction electrons. These electrons 
create a hyperfine magnetic field at the site of the muon, 
which can be considered as a fluctua ting local field with a 
correlation time t c ~ 10~ 13 s in Al |25j. This interaction 
results in a Korringa depolarization rate |2fj, |2jj that has 
an exponential form, and does not depend on the mag- 
netic field. Significant depolarization rates of A > 0.001 
/is -1 have been measured for muons in Cd, Sn, Pb, As, 
Sb, and Bi [2^]. The authors explained the measured 
depolarization rates to be due to Korringa depolariza- 
tion because the A values increase with temperature as 
predicted. 

Jodidio et al. |l2| measured a depolarization rate of 
(0.43±0.34) x 10~ 3 fis- 1 for their Al target at 1.1 T. This 
rate is about 2.5a smaller than the (1.55 ± 0.28) x 10~ 3 
/is -1 observed in this experiment. The difference could 
partly be due to the 2.5 to 5.5% of the muons that stop in 
the gas before the stopping target. The functional form of 
the depolarization in gases is unknown, thus an assump- 
tion that there is no unseen rapid depolarization is made. 
The spin relaxation, with mean lattice-site residence time 
t c , is given approximately by the Kubo-Tomita expres- 
sion |2J|, which reduces to Gaussian (exponential) forms 
for r c — > oo(t c — > 0). 

The difference between Gaussian and exponential ex- 
trapolations of the integral asymmetry measurement, as 
shown in Fig. [3 is 2.4 x 10~ 3 . Data before 1.05 /is are 
not considered because of possible contamination of late 
TDC signals from muons for upstream decay positrons. 
Changing the cutoff between 900 ns and 1100 ns has neg- 
ligible effect on the asymmetry measurement extrapo- 
lated to zero time. Half the difference between the two 
different extrapolations is the correction applied to the 
simulation to data fits, because the simulation was gen- 
erated with a Gaussian form, while in reality the shape is 
most likely a linear combination of a Gaussian and expo- 
nential. An estimate of the extrapolation uncertainty is 



1.005r 



a< - 
| 0.995- 



exponential 
X 2 /NDF = 20.0/20 
Rel. P |1 (0)= 1.0015 ±0.0010 

(1.6±0.3)x 10" f 'ns-' 



0.99 
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Gaussian 
X 2 /NDF = 21.2/20 
Rel. P^O) = 0.9990 ± 0.0006 

o = 53000 ± 5000 ns 
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FIG. 5: (colour online) Extrapolation to zero decay time 
of relative muon polarization, estimated using the decay 
positron integral asymmetry described in the text. The ex- 
trapolation function is fit to data: with an exponential shown 
as a dashed-line, and as a Gaussian shown as a solid-line. 



Other systematic uncertainties were studied by em- 
ploying the fitting technique described in the data analy- 
sis section. In this case the fits are of experimental data 
(or simulation) samples, taken with a systematic param- 
eter set at an exaggerated level, to data (or simulation) 
taken under ideal conditions. The difference measured, or 
sensitivity, expresses the changes in the spectrum shape 
caused by the systematic effect in terms of the changes in 
the Michel parameters. Systematic uncertainties in the 
measurement of P™£, are summarized in Table ITTT1 

The chamber response systematic uncertainty in 
comes from the uncertainty due to the PC and DC re- 
sponse. Six sources of uncertainty, described in the 
following paragraphs, are: time variations in the wire- 
to-wire timing, changes to the space-to-time relations 
(STR) from density variations, chamber foil bulge due 
to changes in differential pressure between the chamber 
gas and cradle helium, the asymmetry of the wire posi- 
tions relative to the foils due to construction, and a dead 
zone due to a recovery time after a muon passes through 
the chambers. 

To estimate the possible time variation in the wire to 
wire time zero (to), calibration runs are taken at the be- 
ginning and end of the run period. To measure the sen- 
sitivity of the muon decay parameters to to variations, a 
calibration scaled by ten times the difference in the be- 
ginning and end of run time calibrations is used. The 
contribution to the chamber response systematic uncer- 
tainty in P££ due to to variations, obtained from fitting 
the exaggerated time shift analysis spectrum to a nomi- 
nal spectrum, is 0.89 x 10~ 3 . A second calibration scaled 
by 5 times the difference in beginning and end of run, 
shows that this systematic uncertainty scales linearly. 

The shape of the chamber cathode foils depends on 
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the differential pressure between the chamber gas and 
the He gas of the cradle holding all of the chambers. The 
differential pressure of the chambers is monitored and 
controlled by the gas system. To measure the level of 
foil bulge, several nominal data runs are taken at dif- 
ferent differential pressures. The average helix fit drift 
time difference, for tracks that go through the center of 
the chambers, relative to the average helix fit drift times 
for tracks that go through the radially distant part of 
the chamber, is sensitive to the foil bulge. Run to run 
monitoring of this fit time difference shows that the vari- 
ation in foil bulge is controlled to better than 50 /um. To 
measure the sensitivity of the muon decay parameters to 
the chamber foil bulge, a simulation generated with the 
STR calibration changed as if the foils were moved out- 
ward by 500 fim is used. The contribution to the cham- 
ber response systematic uncertainty, in P^£, due to time 
variations in the foil bulge, is estimated to be 0.22 x 10~ 3 . 

Surveys of the constructed chambers show that the 
spacing between the wire plane and cathode foils is not 
the same on both sides of the wire. On average the cath- 
odes are found to be shifted by about 150 /im from their 
nominal 2 mm spacing from the wire planes. With the 
solenoid field on, the position of the cathode foils was also 
determined using normal decay positron data. Using the 
bulge calibration data, from runs with different bulges, 
a relationship between a fit time difference to a foil shift 
is obtained. STR files were changed half plane by half 
plane to match the measured shifts. A fit of the muon 
decay distribution, from the shifted cathode foil simula- 
tion to a nominal simulation, results in a contribution to 
the chamber response systematic uncertainty in P^£ of 
0.2 x 10~ 3 . 

History plots of the density variation for each data set 
show that the largest RMS change in density is ±0.7%. 
To estimate the effect of density changes, a simulation 
was generated with the temperature changed by 10% 
(from 300 K to 270 K). A scaled fit of the changed tem- 
perature simulation set to a nominal simulation set is 
used to estimate the addition to the chamber response 
systematic uncertainty in P^£ of 0.17 x 10~ 3 . 

As the muon passes through the upstream half of the 
detector, it deposits large ionization in the chambers. 
The resulting electron avalanches near the wires take 
time to recover as the positive ions drift back to the 
cathode planes. These dead zones produce an upstream- 
downstream efficiency difference for decay positrons. For 
each muon hit the wire efficiency versus distance away 
from the muon hit is estimated as a function of time af- 
ter the muon track. A determination of the muon decay 
parameter sensitivities to the muon dead zone is made by 
introducing an exaggerated dead zone into the simulation 
which removes 11.4 times as many hits as are expected in 
normal data. A conservative estimate of the contribution 
to the chamber response systematic uncertainty in P^£ 
due to the muon dead zone is 0.01 x 10 -3 . 



The systematic uncertainty in P^£ due to spectrome- 
ter alignment comes from how well the chamber transla- 
tions in x, y, z and angle are corrected to match their 
true positions, and the degree to which the misalignment 
between the magnetic field axis to the chamber axis is 
treated. 

The systematic uncertainty in P^£ due to positron in- 
teractions includes four effects: a discrepancy between 
the simulation and data energy loss, hard and intermedi- 
ate interactions, multiple scattering, and backscattering 
from material outside the detector. 

The systematic uncertainty in P^£ due to the momen- 
tum calibration comes from two contributions. One con- 
tribution is due to how well the endpoint energy calibra- 
tion can be determined. The other factor is how well the 
measured magnetic field map represents reality. 

The systematic uncertainty in P^£ due to a differ- 
ence in the upstream and downstream efficiencies is es- 
timated from the mean number of hits used in fitting 
the upstream and downstream tracks. An upstream- 
downstrcam difference is expected from the different an- 
gular distributions, and is reproduced by the simulation. 
The observed discrepancy between data and simulation 
of 0.18 hits/track is a measure of the efficiency difference. 

To estimate the sensitivity of the muon decay pa- 
rameters, due to a difference in upstream versus down- 
stream efficiency, a data set was analyzed with 5% of its 
downstream DC hits thrown away before analysis. This 
upstream-downstream efficiency change, compared to a 
normal analysis of the same data, produces a change 
in P^£ of (1.9 ± 0.9)xl0~ 3 . The lowered efficiency run 
downstream track fits have 1.8 fewer hits/track relative 
to the standard run. An estimate of the systematic uncer- 
tainty in P^£ due to a difference in upstream and down- 
stream efficiency is ergo 0.2 x 10~ 3 . 

The systematic uncertainty in P^£ due to background 
muon contamination comes from muons which are born 
in pion decay in the vicinity of the gas degrader. 

The beam intensity systematic uncertainty in P££ 
comes from changes in the beam rate which affects the 
rate of coincident particles. 

Theoretical uncertainties in the radiative corrections 
introduce a small systematic uncertainty in P^£. 

Several of the systematic uncertainties could vary from 
data set to data set and are denoted by (ave), and 
are considered data set dependent when calculating the 
weighted average value of P^£. For example, the effect of 
positron interactions on upstream and downstream decay 
positrons changes when the mean muon stopping location 
is adjusted; thus the systematic uncertainty in P^£ due 
to positron interactions is set-dependent. 
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TABLE III: Contributions to the systematic uncertainty for 

P£t 

Effect Uncertainty 

Depolarization in fringe field (ave) 0.0034 

Depolarization in stopping material (ave) 0.0012 

Chamber response (ave) 0.0010 

Spectrometer alignment 0.0003 

Positron interactions(ave) 0.0003 

Depolarization in production target 0.0002 

Momentum calibration 0.0002 

Upstream-downstream efficiency 0.0002 

Background muon contamination (ave) 0.0002 

Beam intensity (ave) 0.0002 

Michel parameter rj 0.0001 

Theoretical radiative corrections 0.0001 



V. RESULTS 

The result for presented here uses a data sam- 
ple consisting of 2 x 10 9 events recorded in Fall 2004. 
This data sample includes eight data sets, of which seven 
were used for the extraction of PJJ£. Simulations to fit 
each of the seven data sets used different beam char- 
acterization profiles, derived from beam measurements 
performed after the data collection, which matched dif- 
ferent conditions under which the data were recorded. 
The remaining data set was used to determine the detec- 
tor response using decay positrons from muons stopping 
in the trigger scintillator and the first few chamber planes 
(far upstream), as described in (iqUTi). 

Five sets of data were taken with the beam steered 
nominally. One data set had the muon beam stopping 
with the Bragg peak centered in the target (stop i). Two 
sets, which were separated in time by a few days, were 
taken with the muon Bragg peak shifted to 3/4 of the 
way through the Al stopping target (stop | A, B). One 
set was taken with a muon beam size limiting aperture 
(aperture), and one set was taken with the beam rate 
increased (high rate). 

Two sets of data were collected with the beam dis- 
placed by changing the last bending magnet (B2) field by 
+0.5% from nominal. One of the data sets (B2+0.5%) 
had the muon Bragg peak centered in the stopping tar- 
get, while in the other set (PC5 stop), the muons were 
stopped relatively far upstream in order to increase the 
relative fraction of muons stopping in gas. All of these 
data sets, using different beam characterization profiles 
that matched the different conditions, were used in this 
determination of P^£. 

The spectrum fit results for the parameter P^£ are 
presented in Table llVl At the present stage TWIST can- 
not provide an improved measurement of 77, therefore its 



TABLE IV: Results for P™£. Each fit has 1887 degrees of 
freedom. Statistical and set-dependent systematic uncertain- 
ties are shown. A description of the data sets is in the text. 



Data Set 


P£(, ± stat ± syst 


x 2 


B2+0.5% 


1.0023 ± 0.0015 ± 0.0037 


2007 


PC5 stop 


1.0055 ± 0.0030 ± 0.0038 


1906 


stop | 


1.0015 ± 0.0014 ± 0.0037 


1876 


stop § A 


0.9961 ± 0.0014 ± 0.0037 


1900 


high rate 


0.9997 ± 0.0019 ± 0.0037 


1932 


aperture 


0.9978 ± 0.0018 ± 0.0037 


1896 


stop | B 


1.0009 ± 0.0019 ± 0.0037 


1841 



value is set to the global analysis value of —0.0036 [29j, to 
constrain the other parameters better. The uncertainty 
of ±0.0069 on the accepted value of rj gives an uncertainty 
of ±0.0001 on the final value of P*£. 

The average values of p and 5 from the present fits are 
0.749 and 0.753, respectively. An evaluation of the un- 
certainties in p and 5 has not been performed, but if one 
assumes systematic uncertainties similar to the previous 
TWIST measurements, these values are reasonably con- 
sistent with the published values of p = 0.7508 and 
5 = 0.7496 [13 ■ 

To illustrate the quality of the fit, and how the spec- 
trum fit distinguishes between P^£ and P££S, the contri- 
bution to the fit asymmetry versus momentum for each 
of these terms and from the best fit A{jp) are shown in 
Fig. El Note that the total asymmetry versus momen- 
tum, A(p), is: 

A(p)=At(p)+A cs (p), (8) 

where A% is the asymmetry when the P^S cos 0|(4x 3 — 
3x 2 ) contribution to the positron decay spectrum is zero, 
and A^s is the asymmetry when the P^£ cos #(x 2 — x 3 ) 
contribution is zero. The top panel in Fig. shows 
the best fit asymmetry versus positron momentum, A(p) 
with all of the fiducial cuts applied as a solid line; the 
contribution to the fit from the £ term as a long-dashed 
line; and the contribution to the fit from the £5 term 
as the short-dashed line. The bottom panel shows the 
difference, AA(p) between data and fit. 



VI. CONCLUSION 

The value of P^£ was determined to be 1.0003 ± 
0.0006(stat.) ± 0.0038(syst.). The central value for P*£ 
was calculated as a weighted average using a quadratic 
sum of the statistical and set-dependent uncertainties 
for the weights. The final systematic uncertainty is a 
quadratic sum of the set-independent and the average 
values of the set-dependent systematics. 



Depolarization in fringe field (ave) 
Depolarization in stopping material (ave) 
Chamber response (ave) 
Spectrometer alignment 
Positron interactions(ave) 
Depolarization in production target 
Momentum calibration 
Upstream-downstream efficiency 
Background muon contamination (ave) 
Beam intensity (ave) 
Michel parameter rj 
Theoretical radiative corrections 
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FIG. 6: (colour online) The top panel shows the fit asymme- 
try versus positron momentum, A(p), along with the contri- 
butions to the fit A(p) from the £ and £8 terms. The bottom 
panel shows the difference between the data and fit, AA(p). 
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FIG. 7: (colour online) The top panel shows the manifest 
LRS model 90% confidence limits on £ and mi (gL/gR = 1) 
from measurements of P^£- The bottom panel shows the 
same limits in the general LRS model case. 



The 90% confidence limits on the LRS model parame- 
ters from this measurement are: —0.050 < C < 0.041 and 
m,2 > 360 GeV/c 2 in the manifest case, and —0.061 < 



9r/9l( < 0.061 and gz,/ff.RW2 > 325 GeV/c 2 in the gen- 
eral LRS model. The LRS model limits are shown in Fig. 

HI 

The central value measured is closer to the SM value 
than previous direct measurements, and, hence, in a 
global fit with all other muon decay parameter data [29| 
it pulls those parameters that are sensitive to (Qrr, 
Qlb) closer to the SM value. However, the changes are 
small compared to the uncertainty on these parameters. 
The present result reduces the uncertainty on the direct 
measurement of H by a factor of two; it is also con- 
sistent with the SM and the value obtained indirectly 
0QII3- This is TWIST's first independent measure- 
ment of P5£, and prospects for reducing the main sys- 
tematic uncertainties in P^£ for data taken in the future 
with an improved beam measurement device are excel- 
lent. 
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